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The microscopic global nucleon—nucleus optical model potential (OMP) proposed by Whitehead, Lim, and 
Holt, the WLH potential [1], which was constructed in the framework of many-body perturbation theory with 
state-of-the-art nuclear interactions from chiral effective field theory (EFT), was tested with (p, d) transfer re- 
actions calculated using adiabatic wave approximation. The target nuclei included both stable and unstable 
nuclei and the incident energies reached 200 MeV. The results were compared with experimental data and pre- 
dictions using the phenomenological global optical potential of A. J. Koning and J. P. Delaroche, the KD02 
potential. Overall, we found that the microscopic WLH potential described the (p, d) reaction angular distri- 
butions similarly to the phenomenological KD02 potential; however, the former was slightly better than the 
latter for radioactive targets. On average, the obtained spectroscopic factors (SFs) using both microscopic and 
phenomenological potentials were similar when the incident energies were below approximately 120 MeV. 
However, their difference tended to increase at higher incident energies, which was particularly apparent for the 


doubly magic target nucleus *°Ca. 
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I. INTRODUCTION 


Single-nucleon transfer reactions such as the (d,p) and 
(p, d) reactions are valuable tools for studying the single- 
particle structure of atomic nuclei. In recent years, many 
studies on transfer reactions with radioactive beams have con- 
tributed significantly to our understanding of the evolution of 
nuclear structures and nuclear astrophysics (see, for instance, 
the review papers of Refs. [2—4]). Nuclear reaction theories 
are essential for mediating nuclear structure information and 
measuring nuclear reaction data. In nearly all the theories 
for transfer reactions, for instance, the distorted wave Born 
approximation (DWBA) [5, 6], adiabatic wave approxima- 
tion (ADWA) [7, 8], continuum discretized coupled-channel 
method (CDCC) [9-11], and Faddeev equation-based models 
[12-14], optical model potentials (OMPs) are indispensable 
model inputs. These OMPs affect the theoretical cross sec- 
tions, thus affecting the informal nuclear structure obtained 
from experimental data [15]. Therefore, reliable OMPs are 
crucial for the nuclear reaction calculations. 

Phenomenological global OMPs such as CH89 [16], 
KDO2 [17], Becchetti-Greenless [18], the more recent 
WP [19] for nucleons, and the OMPs cited in refs. [20-25] 
for deuterons or other nuclei are deduced by fitting the elastic 
scattering data of a particle within a wide range of incident en- 
ergies and target masses. Phenomenological global OMPs are 
widely used for nuclear reaction calculations [26-34]. How- 
ever, most experimental data concern stable nuclei. There- 
fore, caution is necessary when applying these phenomeno- 
logical OMPs to reactions involving exotic isotopes. Micro- 
scopic optical potentials, which are based on the more funda- 
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mental principles of nuclear many-body interactions, can be 
more appropriate for describing elastic scattering and trans- 
fer reactions simultaneously for unstable nuclei. Microscopic 
OMPs for nucleons, such as the Bruyères Jeukenne—Lejeune— 
Mahaux model potential (the JLMB potential) [35] and the 
global potential based on the Dirac—Brueckner—Hartree—Fock 
approach developed by R. R. Xu et al. (CTOM [36, 37]), 
and those for deuteron [38], triton [39], a[40], and heavy- 
ion [41], among others, are used with increasing frequency 
in applications to direct nuclear reactions. In addition, mi- 
croscopic OMPs have been found to be more reliable for ex- 
tracting nuclear structure information from transfer reaction 
data [42]. However, many microscopic OMPs are strictly 
semi-microscopic because some parameters constrained by 
nucleon scattering data are still present. 


It has been found that nucleon elastic scattering and trans- 
fer reactions are sensitive to different regions of the OMPs 
[42]. It is possible that phenomenological OMPs that were 
constrained only by elastic scattering data (thus, only a certain 
radial range was well constrained by the experimental data) 
are not the best for theoretical calculations of transfer reac- 
tions. Recently, Whitehead, Lim, and Holt constructed a mi- 
croscopic global nucleon—nucleus optical potential based on 
an analysis of 1800 isotopes in the framework of many-body 
perturbation theory with state-of-the-art nuclear interactions 
from chiral effective field theory (EFT) [1]. An attractive fea- 
ture of the WLH potential is that none of its parameters are fit- 
ted to the nucleon-nucleus scattering data. One might expect 
that being derived fully microscopically, the WLH potential 
might be more suitable for probing nuclear structure informa- 
tion via transfer reactions. Thus, it is necessary to test the 
WLH potential with (p, d) transfer reactions to compare the 
results with experimental data and with the same calculations 
using phenomenological global nucleon—nucleus potentials. 


In this study, tests were conducted using 28 sets of 
(p, d) reaction angular distributions on 15 nuclei, including 
both stable and unstable nuclei, for incident energies below 


200 MeV/nucleon. The global phenomenological nucleon- 
nucleus potential of Koning and Delaroche (KD02 potential) 
[17] was used for comparison. The Johnson—Soper adiabatic 
wave approximation (ADWA) was adapted for the (p, d) and 
(d, p) reaction calculations. The ADWA, which is essentially 
a three-body nuclear reaction model, considers the deuteron 
breakup effect simply but effectively [7]. It is widely used in 
the analysis of (p, d) and (d, p) reactions (see Refs. [43, 44]). 
More importantly, it does not involve the deuteron-target 
OMPs; instead, it considers the deuteron as a p+n system 
and calculates the distorted waves of the (p+n)+target three- 
body system with proton- and neutron-target OMPs. Thus, 
only nucleon-target OMPs are required in the ADWA model, 
which suits our need to estimate the microscopic WLH po- 
tential with (p, d) and (d, p) reactions. 

The remainder of this paper is organized as follows: in Sec. 
II, we briefly introduce the theoretical model and procedures 
of data analysis. An analysis of the transfer reactions is pre- 
sented in Sec. II A. The results for incident energies higher 
than 150 MeV are presented in Sec. III B and for unstable nu- 
clei are given in Sec. II C. A discussion of the spectroscopic 
factors (SFs) is presented in Sec. III D. Finally, we summarize 
our conclusions in Sec. IV. 


Il. METHODOLOGY 


To study the effects of microscopic WLH potential on the 
cross section of the transfer reaction, we investigated the 
(p, d) transfer reactions with 15 target nuclei, which included 
285i, 34,46 Ar, 40Ca, 54Fe, 56,58,60Nj, 907r, 120Sn, 102Ru, 
140Ce, 142,144N dq, and 2°8Pb. The incident energies varied 
between 18 and 200 MeV/nucleon. The ranges of the tar- 
get nuclei and incident energies were mainly limited by the 
upper bounds of the ranges of incident energies of the WLH 
and KD02 potentials, the lower bound of the target mass of 
KD02, and the availability of experimental data. All experi- 
mental data were obtained from the nuclear reaction database 
EXFOR/CSISRS [45] or digitized from original papers, as 
shown in Tables | and 2. 

We adopted the three-body model reaction methodology 
(TBMRM) proposed by Lee et al. to analyze (p, d) reac- 
tions [43, 44, 72]. This methodology uses the Johnson—Soper 
ADWA model for the (p, d) and (d, p) reactions [7], where 
the amplitude of the A(p, d)B reaction is given by [73] 


Tpa = Bm (OSB onplVaplxSy da: (1) 


where S Fpj is the spectroscopic with n, l, and j being the 
node number, angular momentum, and total angular momen- 
tum, respectively, of the single neutron wave function ¢,,; in 
the nucleus A (A = B + n); xpa and xag are the entrance- 
and exit-channel distorted waves, respectively; and V,,,, is the 
neutron—proton interaction that supports the bound state of 
the n-p pair dn,» (the deuteron wave function). 

In this work, Vnp is the Gaussian potential with a depth of 
72.15 MeV and a radius of 1.484 fm, which is taken from Ref. 
[74]. With this potential, only the s-wave was considered in 


TABLE 1. Spectroscopic factors (SFs) obtained with the WLH 
and KD02 potentials for stable isotopes at different energies. The 
reference to the experimental data set used in the extraction is also 
listed. 


Target Ep/MeV Ref SFw iu SFxpo2 
28Si 33.6 [46] 3.353 3.495 
51.93 [47] 3.491 3.727 
185 [48] 1.928 1.419 
“Ca 27.5 [49] 2.570 2.666 
65 [50] 2.568 2.809 
156 [51] 4.364 2.530 
185 [52] 4.482 2.183 
200 [53] 4.114 1.794 
54Fe 51.93 [54] 2.724 3.179 
122.4 [55] 3.187 3.241 
58Ni 24.6 [56] 1.222 1.270 
51.93 [57] 1.113 0.992 
65 [58] 1.251 0.870 
SON; 94 [59] 1.091 0.853 
907% 58 [60] 4.855 5.250 
90 [61] 5.296 6.948 
121.2 [62] 5.772 6.234 
102Ry 26.3 [63] 1.803 1.953 
120S) 18 [64] 2.033 1.903 
26.3 [65] 0.587 0.619 
140Ce 55.05 [66] 2.289 2.029 
142Nd 52 [67] 9.358 8.573 
144Nd 22 [68] 2.147 2.154 
208 Ph 22 [68] 0.703 1.174 
65 [69] 7.018 4.662 

TABLE 2. SFs obtained with the WLH and KD02 potentials for 


unstable isotopes at different energies. The reference to the experi- 
mental dataset used in the extraction is also listed. 


Target Ep/MeV Ref SFwiu SFxpo2 
ST Ar 33 [70] 0.876 0.787 
46 Ay 33 [70] 4.615 5.541 
56Nj 37 [71] 5.380 7.502 


Ọnp- Using the ADWA model, exit-channel distorted waves 
are generated with the following effective “deuteron” (as a 
subsystem composed of neutrons and protons) potential [7, 
8]: 


(ØnplVap [Uns (Ë+ §) + Upe (Ë — §)] lónp) 


(Pnp (T)|Vnp (7) Pnp (T) D 

where Ung and Uppg are the neutron and proton OMPs on 
the target nucleus B evaluated at half of the deuteron inci- 
dent energies (the “E4/2 rule”). In the zero-range version of 
ADWA, the effective deuteron potential becomes 


Uap(R) = 


Uap(R) = Ung (R) + Upp (R). (3) 


Clearly, with the three-body ADWA model, only the neutron 
and proton OMPs are required. This allows the calculation 
of the (p,d) and (d, p) reactions using only one set of sys- 
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Fig. 1. (color online) Comparisons between the two optical model calculations and experimental data of (p, d) reactions on 78Si, “°Ca, **Fe, 
and °SNi at the incident energies, up to 150 MeV, indicated in the figures. The solid and dashed curves are calculated with the WLH and 
KD02 potentials, respectively. The values of the differential cross sections are multiplied by the corresponding SF as labeled. 


tematic OMPs. Zero-range ADWA was found to be satisfac- 
tory for describing the deuteron pickup and stripping reac- 
tions [43, 44]. Therefore, although the finite-range version 
of the ADWA model and more sophisticated yet complicated 
continuum-discretized coupled channel models are available 
for the analysis of (d,p) reactions [8, 73], we adopted the 
zero-range ADWA model in this study. To examine the micro- 
scopic WLH potential and compare its results with the phe- 
nomenological KD02 potential, zero-range ADWA should be 
sufficient. 


The single-particle wave functions were calculated using 
the separation energy prescription with the Woods—Saxon 
form of the single-particle potentials. The depths of these 
potentials were adjusted to reproduce the neutron separation 
energies in the ground state of the target nuclei. The radius 
and diffuseness parameters of these potentials, rg and ao, 
are also important for the nuclear transfer reactions. Their 
empirical values were rọ = 1.25 fm and aọ = 0.65 fm. 
However, these empirical values do not represent the spe- 
cific structure of a single nucleus. A better solution is to 
confine the ro and ao values with a reliable nuclear struc- 
ture theory. The TBMRM constrains the ro and ap values us- 
ing modern Hartree-Fock (HF) calculations [43, 72, 75-83]. 


Using this procedure, the diffuseness parameter was fixed at 
ao = 0.7 fm. The radius parameter rọ was determined by 
requiring the root mean square (rms) radius of the single neu- 
tron wave function, \/(r?), to be related with the rms radius 
of the corresponding single particle orbital from HF calcula- 
tions, \/(r2\Hp, by (r?°) = [A/(A — 1)](r?)up. The factor 
[A/(A — 1)] was used to correct the fixed potential center as- 
sumption used in the HF calculations, where A is the mass 
number of the composite nucleus. All the HF calculations 
performed in this work were based on the SkX interaction 
[84]. We adopted the same procedure as that used in our pre- 
vious studies [42, 81]. All transfer differential cross sections 
were calculated using the TWOFNR code [85]. 


Ill. RESULTS AND DISCUSSIONS 
A. Transfer Reactions on Stable Nuclei at Low Energies 


The reliability of the WLH potential for the (p, d) trans- 
fer reactions was first checked with stable nuclei at incident 
energies Æp < 150 MeV, which fall within the range of in- 
cident energies of the WLH potential. The target nuclei in- 
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Fig. 3. (color online) Same comparison as that in Fig. 1 but on '*°Sn, ‘4°Ce, 4?:!44Nd, and ?°8Pb. 


cluded 28Si, 4°Ca, 54Fe, °8:°Ni, 90Zr, !2°Sn, 192Ru, 140Ce, 
142,144Nd, and 2°8Pb. The results are shown in Figs. 1-3 to- 
gether with the experimental data. For comparison, the results 
calculated using the global phenomenological KD02 potential 
are also presented. All the calculated results were normal- 
ized to the experimental data at the first peaks of the angular 
distributions or at the angles where the maximum measured 
differential cross sections occurred. Neutron SFs were ob- 
tained using this procedure. All experimental data correspond 
to neutron transfer from the ground states of the target nuclei 
to the ground states of the residual nuclei; thus, the results ob- 


tained in this work are single-neutron SFs in the ground states 
of the target nuclei. 


As shown in Fig.1-3, generally, the experimental data can 
be satisfactorily reproduced by the WLH potential at for- 
ward angles as well as by the KDO2 potential. The WLH 
potential was even better than the KDO2 potential in some 
cases, such as '*°Sng. (p, d)'!®Sngs. at 18 MeV. In some in- 
dividual cases, the experimental angular distributions were 
not reproduced satisfactorily using either OMPs, for exam- 
ple, '?°Sn,, (p, d)'!9Sng., at 26.3 MeV. Problems other than 
the OMP may exist in these cases. 
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Fig. 4. (color online) Comparisons between two optical model cal- 
culations and experimental data of (p, d) reactions on 28 Si and *°Ca 
at the incident energies >150 MeV indicated in the figures. The 
solid and dashed curves are calculated with the WLH and KD02 po- 
tentials, respectively. The values of the differential cross sections are 
multiplied by the corresponding SF as labeled. 


B. Extrapolation to Higher Energies 


The WLH potential is expected to work for incident ener- 
gies below 150 MeV. Above this energy range, the theoretical 
uncertainties may become uncontrolled. However, it is in- 
teresting to observe how this operates when extrapolated to 
higher energies. When the incident energies are higher than 
approximately 150 MeV, the randomly generated WLH po- 
tential may generate unphysical positive potential. In these 
cases, we set the number of random pulls to 1500, which al- 
lowed us to obtain the least 1000 negative-valued potentials 
from which we could obtain reasonable averaged potentials. 
In Fig. 4, we show results for (p, d) reactions on ?°Si and 
4°Ca targets at incident energies between 150 and 200 MeV. 
Again, it can be observed that the WLH potential reasonably 
reproduced these higher-energy data and yielded results very 
close to those of the KDO2 potential at forward angles. 


C. Transfer Reactions on Unstable Nuclei 


There is a pressing need for high-quality optical potentials 
for nuclei that are far from stable, which represent a frontier 
in nuclear physics [86-88]. Phenomenological potentials are 
most suitable for reactions where the masses of the target nu- 
clei are within the mass range of the experimental database 
from which the potential parameters are constrained. Be- 
cause most elastic scattering data used to obtain phenomeno- 
logical OMPs were measured on stable targets, the resulting 
OMP parameters may be biased by the limited range of tar- 
get masses in the database. Therefore, caution is always ad- 
vised when extrapolating these potentials to unstable nuclei, 


particularly when they are far from the £ stability line. Mean- 
while, microscopic OMPs, which are derived from effective 
nucleon-nucleon interactions with reliable nuclear structure 
models, are expected to be more sophisticated when applied 
to reactions with unstable nuclei. It is interesting to examine 
how the microscopic WLH potential works on unstable nu- 
clei, and how it compares with the phenomenological KD02 
potential. The results are presented in Fig. 5 for (p, d) re- 
actions on unstable targets °44°Ar and °°Ni. As can be ob- 
served, again, the WLH potential reproduced the experimen- 
tal data reasonably well. In the °°Ni case, the WLH potential 
was even slightly better than the KD02 potential. Unfortu- 
nately, the experimental data for (p, d) reactions on unstable 
nuclei are limited. More measurements of the (p, d) reaction 
data on targets further away from the (-stability line are re- 
quired to examine the applicability of the WLH potential. 


D. Spectroscopic Factors 


The OMPs affect not only the angular distributions but 
also the magnitudes of the (p, d) reactions. Both aspects of 
experimental data are important in nuclear reaction studies. 
Comparisons between the theoretical and experimental an- 
gular distributions determine whether the assumed reaction 
mechanism in the theoretical model and their parameters are 
appropriate, whereas the magnitudes of (p, d) reactions de- 
termine the SFs that provide important information about the 
nuclear single-particle structure. As stated above, the SFs are 
extracted experimentally by matching the theoretical and ex- 
perimental angular distributions at the maximum cross sec- 
tions, where the uncertainties of the experimental data and the 
extracted SFs are at their minimum [70] (see Figs. 1-5). The 
SFs from all reactions analyzed in this study are presented in 
Tables 1 and 2. 

Comparisons between the SFs obtained using the micro- 
scopic WLH potential and the phenomenological KD02 po- 
tential are shown in Fig. 6 for their ratios, SFwuy/SFxpo2, 
at different incident energies. The target nuclei included are 
the doubly-magic isotopes *°Ca and ?°*Pb, stable nuclei ?°Si, 
58,60\j 9977 1205n, 102Ru, 140Ce, and 142:144Nd, and unsta- 
ble nuclei 34-46 Ar and °°Ni. Interestingly, when the incident 
energy was less than approximately 120 MeV, the two sys- 
tematic potentials gave, on average, very similar SFs for all 
types of nuclei. However, when the range of incident ener- 
gies increased to approximately 200 MeV, the ratios showed 
an increasing trend. This is especially apparent for the dou- 
bly magic nucleus *°Ca (there is a lack of experimental data 
for 2°°Pb above 100 MeV). When the results in the tables are 
examined, it can be observed that the changes in SFs of 40Ca 
with the KD02 potential are much smaller over the ranges of 
incident energies than those with the WLH potential. 

It is well known that OMPs with doubly magic nuclei do 
not follow the systematics of OMPs established for other nu- 
clei because of the relatively larger excitation energies of their 
first few excited states [16, 22, 94]. The increased disagree- 
ment between the SFs and the WLH and KD02 potentials 
at higher energies may indicate deficiencies in these poten- 
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Fig. 5. (color online) Comparisons between two optical model calculations and experimental data of (p, d) reactions on unstable nuclei, 
including °“Ar, *°Ar, and °°Ni, at the incident energies indicated in the figures. The solid and dashed curves are calculated with the WLH 
and KD02 potentials, respectively. The values of the differential cross sections are multiplied by the corresponding SF as labeled. 
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Fig. 6. (color online) Ratios of experimental ground-state neutron 
SF values between the WLH and KD02 calculations as a function 
of mass number. The dashed line indicates a perfect agreement be- 
tween data and theory. Open diamonds, open circles, and open tri- 
angles denote SF ratios extracted from (p, d) reactions on doubly 
magic nuclei, stable nuclei(not doubly magic), and unstable nuclei, 
respectively. 


tials. Therefore, we compared the descriptions of proton elas- 
tic scattering from "Ca at higher energies. The results are 
presented in Fig. 7 for the *°Ca target. It can be observed 
that, at lower energies, the two potentials describe the ex- 
perimental data similarly. However, at higher energies, the 
phenomenological KD02 potential appears to be better than 
the microscopic WLH potential. However, neither potential 
reproduced the higher-energy data nor the low-energy data. 
Because all other parameters are the same in the transfer re- 
action calculations, this could be the reason for the differ- 
ences among the SFs obtained at both potentials. Note that 
these incident energies are outside the range of WLH, which 
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is for Ep < 150 MeV. This suggests that caution is needed 
when extrapolating the WLH potentials at higher incident en- 
ergies. Meanwhile, improvement of the systematic potential 
to lighter targets and/or higher-energy regions would be inter- 
esting and useful, especially for doubly magic nuclei. 

IV. SUMMARY 


We verified the performance of the new microscopic opti- 
cal WLH potential based on EFT theory in (p, d) reactions. 
In the present study, we performed zero-range adiabatic cal- 
culations for 15 nuclei, covering a wide range of 18—200 
MeV/nucleon. The phenomenological KDO2 potential pa- 
rameters were used for comparison. This pure microscopic 
optical potential effectively reproduced the angular distribu- 
tions of both stable and unstable nuclei as well as the phe- 
nomenological KD02 potential. Our results suggest that the 
microscopic WLH potential can be used as an advanced ap- 
proach for the prediction of transfer reactions on unstable nu- 
clei and for extracting nuclear structure information of exotic 
nuclei. 

Furthermore, we studied the amplitudes of the transfer 
cross sections with the WLH potential. The SFs extracted us- 
ing the WLH potential were close to those obtained using the 
KD02 potential for all types of nuclei below approximately 
120 MeV. However, when the range of incident energies was 
increased to approximately 200 MeV, the ratios showed a sig- 
nificant increasing trend with an increase in incident energy, 
especially for the doubly magic nucleus *°Ca case. Compar- 
isons between the WLH and KD02 potentials in their descrip- 
tions of proton elastic scattering from *°Ca were performed 
up to 200 MeV. The results showed that some imperfections 
existed at higher energies for both potentials. An improve- 
ment in the WLH potential at higher incident energies for 
doubly magic nuclei is expected. Nevertheless, in general, the 
successful application of the WLH potential in (p, d) transfer 
reactions seems encouraging. 
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Fig. 7. (color online) Comparison of the elastic scattering reaction cross section of the WLH calculation (solid lines), KD calculation (dashed 
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